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A HARDWARE SIMULATION OF A LUNAR MIDCOURSE NAVIGATION 


SYSTEM USING STATISTICAL FILTER THEORY AND 
HAND-HELD SEXTANT OBSERVATIONS 
By Jay V. Christensen and Leonard A. McGee 
Ames Research Center 


SUMMARY 


A simulation investigation is reported that considers the application of 
Kalman statistical filter theory and hand-held-sextant observations to lunar 
midcourse guidance and navigation. The specific system under investigation 
consisted of on-board guidance and navigation computations, a computer control 
and display system, and a hand-held sextant. An on-board digital computer was 
simulated using a ground-based digital computer, and the guidance and naviga- 
tion computations were based on statistical filter theory and linear predic- 
tion. Symbolic abbreviation techniques were implemented in the control of the 
simulated on-board computer and to identify computer data. Observations of a 
simulated celestial scene were taken with a hand-held sextant which had been 
flown on the Gemini XII mission. The observer was a qualified Air Force 
navigator, and observations were taken under static conditions and under 
moving-planet and moving-cab conditions. 

Three sets of eight trajectory runs each were considered for the return, 
moon- to- earth, portion of a lunar mission in which inbound observations were 
processed and inbound velocity corrections were simulated. The resulting 
actual trajectories were compared with the on-board estimated trajectories, 
with the desired reference trajectory, and with theoretical reference data to 
obtain systems guidance and navigation performance data. The actual trajec- 
tories were also compared to acceptable reentry corridor requirements. 

The results further confirm theoretical studies based on assumed-error 
models and indicate that an on-board system that utilizes hand-held-sextant 
observation data, processed by statistical filter techniques, and linear 
prediction, has the potential of providing acceptable guidance and navigation 
performance. No serious anomalies or discontinuities were detected in the use 
of statistical filter processing. Operational results obtained by a symbolic- 
alphabetic abbreviation technique were satisfactory. 


INTRODUCTION 


A number of theoretical studies have shown the potential advantages of 
applying statistical filter theory to lunar midcourse guidance and navigation 
(refs. 1, 2, and 3). These studies were based on processing observational 



data that have an assumed theoretical error distribution. Also, a number of 
studies (refs. 4, 5, and 6) have illustrated the potential of using a hand- 
held sextant for space navigation and guidance. 

To investigate the application of both these concepts, and to further 
investigate actual hardware, the present study considered system guidance and 
navigation performance of a specific hardware-computational system. The 
specific system investigated consisted of a simulated on-board digital com- 
puter, a digital computer control and display panel, and a flight-qualified 
hand-held sextant flown on the Gemini XII mission (refs. 7 and 8). On-board 
digital computations and data processing techniques were used. This investi- 
gation is essentially an extension of the research reported in reference 9, 
which investigated a similar system but used data from a theodolite operated 
manually from a fixed base. 

Performance was evaluated for two environmental conditions : (1) ideal 

static conditions, where there is no motion of the cab or of the planet simu- 
lator; and (2) dynamic conditions, where the observations are taken when the 
cab and planet simulator are moving. Errors for such a system could also 
come from a number of other sources, such as sighting through a window, earth- 
landmark tracking uncertainties, velocity correction maneuver and measurement 
errors, shadow variations on the moon. However, such error sources were not 
considered in this study. 

This study was conducted in three phases. Phase I considered performance 
starting with the last 12 observations, through the last velocity correction 
to perigee. The observation data utilized in this phase were taken by the 
observer in a simulated mission sequence. There were no dynamic motions 
involving the celestial scene or the simulated cab. Phase II included the 
the complete return trajectory and utilized observation data obtained from the 
Phase I investigations. Phase III considered the complete return trajectory 
but utilized observation data taken when the simulated planet was moving rela- 
tive to the star scene and the simulated cab was rotating about all three 
axes. For Phase III, the observation data were not taken in a simulated 
mission sequence as in Phase I. The observation data set used in Phase III 
were obtained over a two-week period under dynamic planet and cab conditions 
typical of worst-case conditions expected during the midcourse phase. This 
technique saved considerable time and provided a practical way of investi- 
gating performance under moving line-of-sight and moving-cab conditions. 

Each study phase consisted of eight data runs in which sextant observa- 
tions were processed, velocity corrections were simulated, and the resulting 
trajectories were computed to vacuum perigee to determine system performance. 
Guidance and navigation performance analysis considered accuracy of the tra- 
jectory state estimation, the resulting velocity correction, the deviation of 
the actual trajectory from the reference trajectory, and the ability of the 
system to meet reentry trajectory state requirements at vacuum perigee. The 
guidance and navigation of the system was compared with theoretical perfor- 
mance of a system which uses identical computations and assumes sextant errors 
having standard deviations that bracket the performance of the actual sextant. 
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SYMBOLS 


A(t e ,ti) 

A-i > a 2 5A3 , A 4 

F 

G 

H 

I 

K 

P 

P' 

P 1 > p 2> p 3 > p 4 

Q 

r 

r 

r 

rms 

r 

r rms 

R 

R' 

A 

R 

^VVP 


prediction matrix (6 x 6) relating state at the end point t e 
to that at ti 

submatrices (3 x 3) of A 

matrix (6 x 6) of partial derivatives of the equation of motion 
with respect to the state variables 

guidance matrix (6 x 6) 

matrix (1 x 3) relating observed sextant angle and vehicle state 

identity matrix with appropriate dimensions 

Kalman weighting matrix (6 x l) 

covariance matrix (6 x 6) of estimation error 

new value of P after a sextant observation 

submatrices (3 x 3) of P 

covariance matrix (1 x 1) of the observational errors 
actual position deviation from reference ( 3 x 1 matrix) 
error in estimate of r (3 x 1 matrix) 

rms value of the error in the estimate of r (3 x l matrix) 
estimate of r (3 x 1 matrix) 
rms value of the estimate of r 

covariance matrix (6 x 6) of the deviations between the actual 
and reference states 

new value of R (6 x 6) after a velocity correction 
increment 

Gaussian noise with zero mean 

transition matrix (6 x 6) relating state at t, to state at 
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Notation Conventions 


T 

( ) 

( r 1 

E[ ] 

trace [ ] 

Av rms 

AV 

AV G 

|4V g | 

X 

X 

X 

X,Y,Z 

Y V 7 
A m> ’iri^m 

X s , Y s , Z s 

y 

y 


transpose of matrix ( ) 
inverse of matrix ( ) 
expected value of [ ] 

sum of the principal diagonal elements of [ ] 
rms value of the indicated velocity correction 
covariance matrix (3 x 3) of the indicated velocity correction 
velocity-vector increment (3xl matrix) to be gained 
magnitude of AVq 

six-vector (6 x l) deviation of the actual position and 
velocity from a reference trajectory (x + x) 

estimate (6 x l) of x 

error (6 x 1) in the estimate of x 

Cartesian coordinates of vehicle's position 

Cartesian coordinates of moon's position 

Cartesian coordinates of sun's position 

actual sextant space angle 

estimate of y 


e 

k 

t 

rms 


Subscripts and Matrix Arguments 

end point 

at the kth observation 
general time arguments 
end-point time 

time of the kth observations 
root-mean- square value 
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SIMULATION DESCRIPTION 


Specific hardware and an on-board computer simulation were used to 
investigate on-board techniques for navigation of manned spacecraft. The 
system evaluated in this investigation is illustrated in figure 1. 


Hand-Held Sextant 

The sighting instrument (fig. 2) was a hand-held sextant which had been 
flown on the Gemini XII mission. The general characteristics of the sextant 
follow : 

Size: 7-1/4 x 6-1/16 in. 

Weight: 6-1/2 oz 

Magnification: 8x 

Measurement range: 76° 

Field of view: 7° 

Reference 7 describes the sextant design and space-flight rating. 


Computer Control and Display 

The digital computer control and display panel concept was based on 
symbolic abbreviations. The display panel consisted of the following sub- 
assemblies as numbered in figure 3: (1) electroluminescent computer control 

mode lamps, (2) incandescent status and warning lights, (3) time interrogate 
pushbutton, (4) real-time clock display and control, (5) electroluminescent 
symbolic and numeric display panels, (6) single function, priority interrupt 
pushbutton panel, (7) thumbwheel input panel, and (8) keyboard input panel. 

Details regarding the computer control and display panel concept, opera- 
tion, data format requirements, symbolic code structure, computer memory 
requirements, and functional descriptions of the required programs are 
discussed in reference 9. 


On-Board Digital Guidance and Navigation Computations 

The on-board digital computer was simulated with a general-purpose, 
ground-based computer. The mathematical formulas and techniques for on-board 
guidance and navigation computations were essentially those used by Smith, 
McLean, Schmidt, and McGee (refs. 1, 2, and 3), and are discussed further in 
the section on digital computations. The integration technique, reported in 
reference 10, was a Stormer-Cowell integration with a starter that builds up 
the table of differences backwards in time from the point of integration. 

This particular technique was chosen to reduce integration error as much as 
possible. Knowledge of the position of the sun and moon was required, and a 
Tchebycheff polynomial in time was fitted to each Cartesian coordinate of the 
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sun-moon data over a 17-day time period. This time span was chosen to allow 
for the 7-day round-trip, launch delays, and stay-time on the moon. The 
polynomial coefficients were computed from Naval Observatory data. 

For the on-board computer simulation, two sets of nonlinear equations of 
motion were computed and integrated: the reference trajectory set, which 

defined a nominal trajectory around which the equations of motion were 
linearized for prediction and guidance; and the estimated trajectory set, 
which was computed to obtain the best estimate of the position and velocity 
of the vehicle. The equations of motion included the effect of the sun, the 
moon, and the second harmonic of the earth's gravitational potential; these 
are summarized in appendix A. The coordinate system used was a nonrotating, 
right-hand orthogonal Cartesian frame with its origin at the center of the 
earth, as illustrated in appendix A. 


SIMULATION CONFIGURATION AND TECHNIQUES 


The simulation configuration is illustrated in figure 4. The computer 
control and display panel was mounted in a moving-base cab simulator, and 
observations of the simulated scene were tkaen by the observer under dynamic 
and static conditions. 


Digital Computer Simulation 

An SDS 920 computer was used to simulate the on-board computer and to 
perform the necessary calculations, calibrations, and data processing associ- 
ated with the celestial scene and the computation of the research data. The 
computer system configuration is shown in figure 5. The simulation was a 
complete digital simulation. The digital central processor was a 24-bit, 
general-purpose machine with 12,288 words of random accessible storage. 

FORTRAN and symbolic programming techniques were used. The advantages of 
FORTRAN in meeting the complexity and scaling requirements of the guidance 
and navigation equations outweighed the problems of interfacing the FORTRAN 
with the machine language (symbolic) portions of the program. The real-time 
hardware systems, such as the computer control and display programs, real-time 
mission clock control and display, and celestial scene encoder interrogation, 
could only be controlled by symbolic machine language programs. They were 
written in symbolic language and called through interrupts or as FORTRAN sub- 
routines. Storage of 24,000 words was required for the simulation; of this 
amount, 12,000 were required for the on-board system simulation. Since the 
total program exceeded core storage, it was necessary to fractionalize the 
program into links on the magnetic tape and call in portions as appropriate. 
Double-precision, floating-point programs having 11 significant decimal digits 
were used for all guidance and navigation computations. More details regard- 
ing the digital computer simulation techniques used are discussed in 
references 9 and 11. 
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Digital Computations 


The digital computations used in this study were a modified version of 
those used for the midcourse navigation studies reported in reference 1. The 
reference mission was the return leg of a circumlunar trajectory. Imperfect 
injection conditions for the space vehicle at the perilune point were assumed, 
and the function of the navigation and guidance system was to cause the space 
vehicle to arrive at the reference end point (vacuum perigee) with a minimum 
of time and position error. The navigation system processed sextant observa- 
tions of the angle between a known landmark and a known star at preselected 
times. From the resulting data, the system generated an estimate of the vehi- 
cle state (position and velocity) by means of a Kalman statistical filter. 

On the basis of this estimate, the guidance system executed exactly the 
estimated corrections to the vehicle velocity that were required in order for 
the vehicle to arrive at the reference end point (vacuum perigee) . 

Twenty-seven, second-order ordinary differential equations, comprising 
three trajectories and six sets of three perturbation equations, were inte- 
grated simultaneously. A numerical integration subroutine described in 
reference 10 was used in these computations. 

For this simulation study, three trajectories were computed: (1) a 

previously determined reference trajectory, (2) an estimated trajectory 
giving the current estimate of the vehicle state, and (3) an actual trajec- 
tory giving the true state of the space vehicle. The reference and estimated 
trajectories were processed in the on-board computer simulation. The actual 
trajectory would not be available on board an actual spacecraft, so it was 
considered to have been computed external to the spacecraft. Comparison of 
statistical deviations between the states of these three trajectories 
permitted evaluations of the performance of the system with various sextant 
errors . 

Figure 6 is a block diagram illustrating the computational flow: the 

bracketed numbers in the following discussion identify the steps in the 
digital computation process, and correspond with numbered boxes in the figure. 

In this simulation investigation, when a sextant observation made at 
time t^ is to be processed, numerical integration of the equations of 
motion is initiated simultaneously [1], for the reference, actual, and esti- 
mated trajectories. The equations of motion and the coordinate system used 
are described in appendix A. Initial conditions for the actual and reference 
trajectories are merely the states that existed at the time of the last obser- 
vation, whereas the estimated trajectory uses the new estimated state derived 
from processing the last sextant observation. Simultaneously with the inte- 
gration of the equation of motion by the numerical technique of reference 10, 
the sixfold perturbation equations [2] are also integrated numerically from 
the time of the last observation %_ 1 . This process continues until the com- 
puter time equals the time of the observation tp. The integration is then 
stopped, and the transition matrix < Kt] < ,tj < _ 1 ) is calculated [3] as described 
in appendix B. 


7 



A discussion of the mathematics used in the computations described in the 
remainder of this section is contained in appendix C. 

The covariance matrix of the estimation error P (t^) and the covariance 
matrix of the deviations between the actual and the reference states R(tp) 
are updated [4] through linear prediction. 

Next, the matrix H, relating the observed sextant angle to the vehicle 
state, is computed [5] so that the Kalman weighting matrix Kft^) may be 
obtained. The Kalman weighting matrix is computed [6] for this application 
in the form given in reference 1. With K(t^) known, the value of the covari- 
ance matrix of estimation error after a sextant observation P'ft^) is then 
computed [7]. The delay unit [8] represented the storage of PCt^) until the 
time of the next observation when it becomes P(t^_ 1 ). 

The next process is to compute the new estimated state through a series 
of steps. In this system, a sextant observation consists of measuring the 
angle between a chosen star and a chosen landmark. First, the estimated sex- 
tant angle is computed [9] from geometrical equations using the estimated 
position; then the actual sextant angle is computed [10] using the actual 
position. The estimated sextant angle is next subtracted from the actual 
sextant angle, and observation errors are added, which are based on the 
covariance matrix of the observational errors, Q, or on errors obtained by 
taking observations of a simulated celestial scene with an actual sextant. 
Finally, the result [11] multiplied by Kft^) [12] is added to the estimated 
state [13] to produce a new estimated state. 

Next, the indicated velocity correction is computed by means of a 
prediction matrix Aft^t^), which predicts the end-point state from the 
present state [14]. Submatrices from Aftg,^) are used to compute a matrix 
G which is needed in computing the indicated velocity correction. The indi- 
cated velocity correction, AVq, is computed by multiplying the difference 
between the reference state and the estimated state, x [15] by G [16]. The 
statistics of the indicated velocity correction are determined from a 
covariance matrix AV [17] . 

This concludes the computation cycle at each observation, and the 
computer is now ready to process the next observation unless a velocity correc- 
tion is to be executed. In this case, the correction is made after a time 
delay simulating the time required to orient the vehicle along the desired 
thrust vector. Trajectory determination processes [1], [2], and [3] are con- 
tinued until the computer time equals the time at which the velocity correction 
is to be made. Processes [4], [14], [15], and [16] are repeated, and computa- 
tional control is passed to a separate program section for the performance of 
the velocity correction maneuver. 

The flow lines that represent the velocity correction process are shown 
as dashed lines in figure 6 to emphasize the occasional nature of the veloc- 
ity corrections. Also, these corrections are assumed to occur instantane- 
ously so that no position or time change need be accounted for during the 
maneuver. The velocity correction results in a change in the subsequent actual 
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and the estimated states. To complete the velocity correction cycle, the R 
matrix, which is the covariance matrix of the errors between the actual and 
the estimated states, must be corrected to reflect the effects of the velocity 
correction maneuver [18] . The computation then proceeds to process the next 
observation . 


Celestial Scene Simulation 

Figure 7 is a photograph of a star simulator, which is a collimated 
point light source fixed with respect to the moon simulator base and the 
spacecraft cab simulator base. The moon simulator (fig- 8) projected a colli- 
mated photograph of the moon which subtended an angle of 1/2° from limb to 
limb. The collimated image of the moon was either repositioned between obser- 
vations or driven at a constant angular rate with respect to the fixed-star 
source. The moon simulator interfaced with the computer through a natural 
binary digital encoder that provided an angle readout resolution of 1.25 sec- 
onds of arc. Figure 9 is a view of the simulated celestial scene from the 
observer's station in the cab. In this study, angular travel was restricted 
under dynamic conditions to less than 7 minutes of arc to enhance accurate 
angle generation. 

A set of 96 observations was taken of the moon and star simulator using 
the hand-held sextant under ideal static conditions (Phase I data) . The 
resulting sextant error set had a mean of 2.5 seconds of arc and a standard 
deviation about the mean of 5.9 seconds of arc. A comparison of these 
results with those of reference 6 shows good correlation (5.9 sec error as 
compared with 4 and 5 sec) . One can conclude then that very little random 
error was generated by the simulation of the celestial scene. The 2.5 sec- 
onds mean error may be assumed simulation calibration error or instrument 
calibration error. In both cases (random and mean error) it was felt that 
this magnitude of the simulation error was not large enough to compromise the 
results of the study. 


Spacecraft Cab Simulation 

The spacecraft cab simulator is shown in figure 10. The computer 
control and display panel was installed in the cab so the observer could con- 
trol the simulated on-board computer (fig. 11). The hand-held sextant was 
suspended by a cable from the top of the observation window to reduce the 
fatigue effects of holding the instrument in a gravity field. Particular care 
was taken in this installation to reduce the possibility of adversely stress- 
ing the instrument. When observations were taken, the background lighting of 
the scene which was required to take the photographs of figures 7-11, was 
eliminated. The cab was mounted on an air-bearing to reduce friction and 
rotated about all three axes to simulate a limit-cycle condition for the 
dynamic phase of the simulation. An analog computer was used to generate 
appropriate signals for the servo-drive system that rotated the cab. 
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INVESTIGATION PROCEDURES 


A representative moon-to-earth return mission and trajectory was chosen. 
The selected observation schedule, initial conditions, and system errors were 
typical of a moon-to-earth return guidance and navigation problem; however, 
no attempt was made to study or optimize these factors. The reference tra- 
jectory was the return portion of a free return circumlunar trajectory of a 
spacecraft launched from Cape Kennedy. Closest approach to the moon (peri- 
lune) was 185.3 km (100.0 nautical miles) from the surface. Vacuum perigee 
was 60.0 km above sea level at latitude 24.2° N, and longitude 80.98° W, 
which is approximately over Havana, Cuba. 

All mission times were from trans lunar injection. The spacecraft 
injection at perilune takes place at 70.68 hours. Initial errors were iden- 
tified by adding initial condition errors to the reference trajectory at the 
perilune mission time. A random number generator computed these initial con- 
dition errors on the basis of a normal error distribution with a standard 
deviation of 1 km and 0.001 km/sec in each Cartesian coordinate. The 
specific errors that resulted are given in table I. 


Observation Procedures 

Two types of observation sets were taken. One set was taken under 
static conditions where, prior to each observation, the moon simulator was 
positioned to a new angle unknown to the observer; motion of the celestial 
scene or the cab was not introduced. The other observation set was taken 
under dynamic conditions where cab motion and planet motion were introduced 
during the observation. During these dynamic observation data runs, the line- 
of-sight to the moon simulator was moving with respect to the fixed star 
simulator at 5 seconds of arc per second, which is considered a maximum rate 
for the midcourse phase. 

All observations were taken with respect to the center of the moon 
crater Tycho (fig. 12). All observations and velocity corrections were based 
on fixed times as summarized in table II and illustrated in figure 13. An 
observation was taken by superimposing the star source and the center of the 
moon crater as seen through the sextant. In all cases, observations were 
taken with respect to the moon simulator even when the schedule called for an 
observation of the earth. When the observer was satisfied with the observa- 
tion (i.e., when the targets were superimposed in the sextant sighting task), 
he actuated a mission time-interrogate pushbutton mounted on the sextant after 
which the sextant angle was not modified. The activation of the mission time- 
interrogate pushbutton initiated a computer program that stored the mission 
time of the observation to 0.01 second, interrogated the moon simulator 
digital encoder to determine the moon simulator angular position, and then 
computed the simulated celestial scene angle that existed at the time of the 
observation. The observer read the observed angle from the index readings on 
the sextant and entered it into the computer through the computer- control and 
display panel. This sextant observation angle was then compared with the moon 
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simulator angular position and the difference was taken as the error in the 
observation. When an earth observation was scheduled, the on-board computa- 
tions processed the observations of the moon simulator as if the observations 
had been taken of the earth. 

All three trajectories (estimated, reference, and actual) were then 
updated to the time of the observation, and the "actual" angle that should 
have been observed at that point from the actual trajectory was computed for 
the particular star and planet. The observation error was then algebraically 
added to the "actual" angle, and the resulting "observed" angle was then 
processed using the statistical filter theory. This simulation technique is 
discussed further in reference 9. 


Velocity Correction Procedures 

The velocity correction maneuver was simulated. The velocity correction 
magnitude and direction, as estimated by the on-board guidance and navigation 
computations, were used in the simulation. Errors in cutoff, alinement, and 
measurement were not introduced. 


Study Phases 

The three study phases are summarized as follows: 



Phase I 

Phase II 

Phase III 

Mission segment 

122.0 hr 

70.68 hr 

70.68 hr 

Number of observations 

to perigee 
12 

to perigee 
39 

to perigee 
39 

Velocity corrections 

Last inbound 

All 3 inbound 

All 3 inbound 

Conditions under which 

Static 

Static 

Planet and 

the observations 


(processed 

cab motion 

were taken 


Phase I data) 


Number of data runs 

8 

8 

8 


In Phase I, the observations were taken, and entered, and all guidance 
and navigation computations were carried out on-line, in a real-time sequence 
typical of that phase of the mission. Phase I was restricted to the observa- 
tion sequence immediately preceding the final velocity correction. In 
Phases II and III, observation error sets were processed off-line to allow an 
investigation of a number of complete return trajectories. 

Phase I .- The Phase I study processed each of 12 observations as they 
were taken for the last velocity correction sequence, starting with the obser- 
vation scheduled for 123.0 hours. The trajectory was divided into two parts: 

(1) the generation of the initial conditions at the 122.0 hour point, and 

(2) the study sequence covering the last 12 observations and the final return 
velocity correction. The initial conditions at the 122.0 hour point were 
obtained by starting the inbound trajectory at the moon-to-earth injection 
point (70.68 hr) with injection errors as previously described. A random 
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number generator was then used to generate instrument sighting errors based 
on a gaussian distribution with 10 seconds of arc standard deviation. The 
first 27 inbound observations were computed and processed, and both velocity- 
corrections were simulated through a mission time of 122.0 hours, according 
to the observation and velocity correction schedule outlined in table II and 
illustrated in figure 13. The resulting errors in position and velocity 
obtained at the 122.0 hour point, shown in table III, constitute the initial 
condition errors used in all the subsequent Phase I data runs. 

The real-time mission clock was initialized appropriately so that all 
observations were taken in a simulated real-time mission situation. Prelim- 
inary information that the observer had to enter prior to each observation 
included: observation number identification, star identification, and planet 

identification (moon or earth) . After each observation of the celestial 
scene (as described previously), the actual and on-board trajectories were 
updated to the time of the observation. On the basis of these observations, 
the estimated velocity corrections were computed and then simulated. After 
the velocity corrections were executed, the resulting actual trajectories 
were integrated to the time of the reference trajectory vacuum perigee and to 
the actual vacuum perigee. Eight missions were completed in this manner. 

Phase II . - The Phase II study used the actual observation set previously 
obtained from Phase I in a manner that would permit system performance anal- 
yses covering the entire inbound trajectory. The trajectory conditions were 
initialized at the moon-to-earth injection point identical to Phase I 
(table I) . The 96 observation errors obtained in the Phase I data runs were 
then used by sequential cycling through the set to obtain all the observation 
errors required. These observations were processed, and all scheduled veloc- 
ity corrections were simulated at the times given in table II. After the 
final velocity correction, the resulting trajectories were integrated to the 
time of the reference vacuum perigee and to the actual perigee time. Eight 
missions from injection to perigee were executed in this manner. 

Phase III .- Phase III used observation data taken under dynamic condi- 
tions. The observer took observations while the moon simulator was moving 
with respect to the star background at 5 seconds of arc per second, consid- 
ered a maximum rate for the midcourse phase. At the same time, the cab 
(fig. 10) was driven with a limit cycle of ±2° in all three axes at 0.08 Hz. 
Sixty-seven observations were taken off-line over a two-week period to reduce 
effects of fatigue. The resulting observation errors were used sequentially, 
as described for Phase II. After the final velocity correction, each result- 
ing trajectory was computed to the time of the reference vacuum perigee and 
to the actual perigee time. Eight moon-to-earth return missions from 
injection to perigee were executed in this manner. 


Theoretical Reference Data 

To compute theoretical reference data to which the experimental data can 
be compared, variances of the injection errors are assumed known and normally 
distributed with zero mean. This statistical information in the form of 
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covariance matrices P and R allows the computation of the statistics of 
the entire ensemble of trajectories from one time to the next time of 
interest. The accuracy with which these covariance matrices represent the 
entire ensemble of trajectories depends upon the validity of the linear per- 
turbation equations that are used to compute the transition matrix $ ; and 
the validity of the linear perturbation equations used is discussed in 
reference 1 . 


The 6x6 covariance matrix P (see box 4 of fig. 6, and appendix C) 
represents statistics of the error between the actual trajectory and the 
estimated trajectory, and is a fundamental part of the Kalman filter used in 
processing observations to obtain a better estimate of the position and 
velocity. As such, it is affected by the processing of each observation and 
depends on the assumed statistical performance, which in turn is based on the 
covariance matrix Q of sextant statistics, contained in the covariance 
matrix Q. The P matrix is partitioned as follows: 


P = 


Pi 

?3 


Pz 

Pk 


( 1 ) 


and the rms magnitude of the position uncertainty of the estimated trajec- 
tory with respect to the actual trajectory ? rms is computed as follows: 


r rms 


[trace Pi] 1 ^ 2 


( 2 ) 


The 6x6 covariance matrix R (see box 4 of fig. 6, and appendix C) , 
which allows the computation of the statistics of the entire ensemble of 
actual trajectories with respect to the reference trajectory, is indirectly 
affected by the observation through the velocity correction. The R matrix 
is partitioned as follows: 


R = 


R 1 

R 3 


R 2 

r 4 _ 


(3) 


and the rms magnitude of the position uncertainty of the actual trajectory 
with respect to the reference trajectory r rms is computed from the 
following : 


r rms 


[trace R p ] 


1/2 


(4) 


The rms magnitude of the indicated velocity correction Av rms at any 
particular time along the trajectory depends on both P and R, and is com- 
puted from the 3 x 3 Av matrix (see box 17 of fig. 6 and appendix C) . The 
value of Av rms is used to evaluate the effect of sextant observation errors 
on the velocity correction, and is computed as follows: 
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Av rms = [ trace Av] 1/2 (5) 

Values of r rms , r rms , and Av rms were computed for sextants having 

errors with a standard deviation of 1, 10, and 50 seconds of arc. These 
values were computed with appropriate initial conditions for each of the 
three phases and were used as base line data. The experimental results are 
compared to this theoretical data in the discussion of the results . The 
standard deviation of the errors in the instrument is indicated in seconds of 
arc; for example, f rms (50) refers to data obtained assuming a sextant having 
errors with a standard deviation of 50 seconds of arc. 


RESULTS 


System guidance and navigation performance was evaluated on the basis 
of the following: (1) trajectory estimation accuracy by comparing the esti- 

mated trajectories with the actual trajectories; (2) velocity correction 
control by comparing the velocity correction magnitudes with theoretical 
reference data; (3) trajectory control by comparing the resulting actual tra- 
jectories with the on-board reference trajectory through actual trajectory 
vacuum perigee points, and by comparing data at the actual trajectory vacuum 
perigee with allowable limits. 

Each research phase consisted of eight complete data runs. For 
graphical presentation the experimental results were reduced statistically to 
produce a composite trajectory envelope representing the standard deviation 
about the mean. This envelope, which is plotted for comparison with the 
theoretical reference data, consists of a "system performance 1 sigma upper 
boundary," and a "system performance 1 sigma lower boundary." The upper 
boundary curves were computed by adding the 1 sigma value computed for the 
eight data runs to the computed mean of the eight data runs, and the lower 
boundary curves were computed by subtracting the 1 sigma value from the mean. 

In all the data runs using the experimental data, the Kalman filter 
equations assumed a standard deviation of the sextant error of 10 seconds of 
arc. The actual instrument data obtained were slightly better than those 
modeled in the Kalman filter. Consequently, a slightly mismatched situation 
resulted. Also, a specific, but typical, member of the ensemble was chosen 
to initialize the experimental data; this resulted in slightly degraded ini- 
tial conditions for the experimental data as compared to the statistical 
conditions used in initializing the theoretical reference data. 


Trajectory Estimation Accuracy 

Trajectory estimation accuracy (1 sigma envelope) obtained in Phase I is 
illustrated in figure 14. The accuracy of the system performance estimation 
was generally equal to or better than that theoretically predicted for a sys- 
tem using an instrument with a standard deviation error of 10 seconds of arc. 
This can be seen by comparing the performance envelope with the theoretical 
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reference curve r rms (10). The system estimate after the last observation had 
an accuracy that an interpolation of the reference data indicates would result 
from processing observation data having 7-8 seconds-of-arc error. This tra- 
jectory estimation performance is considered acceptable because the resulting 
actual trajectories provide a safe reentry. 

As previously mentioned, the initial conditions were obtained assuming a 
standard deviation error in the instrument of 10 seconds of arc and a statis- 
tical estimate of the instrument in the Kalman filter of 10 seconds of arc. 
Subsequent to the initialization point (fig. 14, point A), theoretical data 
were obtained for three reference cases: 1 second of arc, 10 seconds of arc, 

and 50 seconds of arc. For the 1 second case t rms (l), the subsequent data 
processing was sensitive to the more accurate measurements, as indicated in 
figure 14 by an immediate improvement in position estimation accuracy. 

For r rms (10) , as would be expected, there was no significant change in 
the estimation accuracy. When the standard deviation of the error in the 
instrument was changed to 50 seconds of arc r rms (50) a significant advantage 
of the Kalman filter was illustrated] the filter gives little weight to these 
degraded observations since the uncertainty of the estimated state (due to the 
higher initialization accuracy) is low compared to the new measurement 
uncertainty . 

Phase II covered the complete mission, using the data of Phase I as 
previously described. As illustrated in figure 15, the system position esti- 
mation accuracy for this phase was very similar to that obtained for Phase I, 
and is considered quite adequate. An interpolation of the reference data 
indicated that performance of the system was equivalent to that indicated by 
the theoretical reference data for an instrument with about 5-6 seconds of 
error (standard deviation). In this phase, the theoretical reference data 
for had time to settle out and illustrate that performance of this 

type of instrument is unacceptable. 

The Phase III data for the dynamic situation, where moving line-of-sight 
planet motion and cab-limit cycle motion are introduced, are illustrated in 
figure 16, which compares system estimation performance data through the 
complete in-bound trajectory with theoretical reference curves. The system 
estimation performance settled out around 20 kilometers of position estima- 
tion error. An interpolation of the theoretical reference curves indicated 
that this level of performance would result from processing observation data 
having 11-seconds-of-arc error, and thus would provide a safe, fixed-time-of- 
arrival reentry, which is consistent with theoretical predictions. The 
results are also consistent with the statistics of the observation error set 
obtained for this phase which, for the 67 observations that were taken, had a 
mean of 4.7 seconds of arc and a standard deviation of 8.6 seconds of arc 
giving an rms value of 9.8 seconds of arc. A comparison of statistical 
errors of the observation sets of Phases I and III indicates that the error 
increase due to the dynamics of Phase III is about 2.2 seconds of mean error 
(4.7 - 2.5), and that the standard deviation error increased by 2.7 seconds 
(8.6 - 5.9) . 
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Velocity Correction 


Table IV compares the velocity correction magnitudes obtained for all 
three phases with theoretical reference data. The comparison of velocity 
correction magnitudes of Phase I with theoretical reference data indicates 
satisfactory performance. There was a trend in Phase I for the velocity 
correction magnitude to be greater for the more accurate sextant models, as 
illustrated by a comparison of Av rms (50) with Av rms (l). This trend was an 
indirect result of the Kalman computation technique. An improvement in the 
knowledge of estimated position and velocity results in an increase in the 
magnitude of the velocity correction. When a series of corrections was made, 
as in Phases II and III, the sum of the velocity corrections showed that the 
more accurate instrument not only provided better actual position control 
(discussed below) but also minimized the expenditure of fuel, as one would 
expect . 

The total velocity correction magnitudes obtained in Phases II and III 
(line 5, table IV) are consistent with the theoretical reference data. The 
magnitudes of Phase III, for example, are slightly larger than those obtained 
for a system having a theoretical 10-second standard deviation of the error 
in the instrument Av rms (10), which is consistent with the position estimation 
performance discussed in the previous section. 


Trajectory Control 

Errors between the actual trajectory and the desired reference trajectory 
are shown for Phases I, II, and III, in figures 17, 18, and 19, respectively. 
Performance for all three study phases shows acceptable guidance and naviga- 
tion. The increase in the upper error envelope boundary as the actual trajec- 
tory approaches the perigee point is to be expected because of the propagation 
of the uncertainties in the velocity correction. This upper error boundary 
still represents adequate systems performance; the actual trajectory errors 
from the reference trajectory result in a safe reentry state, as discussed in 
the following paragraph. Again, a comparison with theoretical data shows the 
resulting performance to be consistently better than that expected for a 
10-second-of-error instrument. 

Inasmuch as the altitude deviation at the actual perigee point is 
particularly critical and the errors in downrange and crossrange are also of 
interest in evaluating fixed time of arrival performance, the actual trajec- 
tory state is evaluated specifically at this point. The data are summarized 
in table V. Performance in position control at the actual perigee point was 
acceptable in all phases, as indicated by a comparison with the reference 
data and the allowable reentry corridor. The allowable reentry corridor data 
for this figure were obtained from references 12 and 13. Most of the position 
error at the actual perigee position is in the downrange direction and repre- 
sents an error in fixed-time-of-arrival guidance of but a few seconds (the 
spacecraft velocity at this point in the downrange direction is approximately 
11 km/sec) . 
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Observer Comments on the Use of the Hand-Held Sextant 


The observer was a qualified Air Force navigator (Capt.) with 
considerable experience with this type of instrument. He felt the hand-held 
sextant was operationally satisfactory for sightings of the type required in 
this experiment. The moon crater Tycho, used for the landmark reference, was 
easy to -locate but there was some difficulty experienced by the observer in 
consistently determining the center accurately. He felt that actual sightings 
on the moon would be more difficult, mainly due to shadows of the craters 
shifting each day. 

There were no significant problems during the static observations. The 
observer did find it desirable to introduce a small amount of instrument rota- 
tion motion during the observation. This motion produced a slow frequency 
motion of the star relative to the landmark with a magnitude as illustrated 
between points A and B of figure 12, and provided a more precise method of 
superimposing the star in the center of the landmark. During the dynamic 
measurements, the line-of-sight rates seemed very slow (5 seconds of arc per 
second) . He felt that taking an observation was not any more difficult due to 
the moving line-of-sight variable than taking the observations under static 
conditions. However, the effects of the moving cab were significant, as this 
motion made it difficult to keep the star and the landmark superimposed. The 
motion of the cab coupled into the sextant and resulted in the same kind of 
star- landmark motion that was introduced by the observer in the static phase. 
However, this effect was far greater than desired. The observer attempted to 
reduce it by rotating the hand-held sextant in a direction opposite to the 
cab motion, but the resulting motion of the star and the landmark could not be 
reduced to less than that indicated by points C and D of figure 12. It is 
felt that this motion was the main source of error increase in the dynamic 
Phase III observations over those obtained from the static observations of 
Phases I and II. 

The guidance and navigation performance obtained in this study were 
adequate, even though the cab motion used (2° at 0.08 Hz in all axes) was con- 
siderably in excess of that expected of an actual control system, and the 
coupling effects of this on the hand-held sextant were greater than desired. 

It is expected that under gravity-free conditions of space and isolation of 
the sextant from physical attachment to the shell of the vehicle, the error 
due to cab motion will actually be less than that experienced in this 
experiment . 


Symbolic Computer Control and Display 

The symbolic abbreviation computer control and display concept used in 
this study was identical to that reported in reference 9, except the panel 
configuration was modified for greater compatibility with cab installation 
requirements. Very satisfactory operational results were obtained through 
application of the symbolic-alphabetic abbreviation technique for entering 
and interpreting computer data. No errors were generated during input that 
were not detected and corrected by the observer, and there were no detectable 
display interpretation errors. The observer became familiar very rapidly 
with the abbreviations and the sequencing. 
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CONCLUSIONS 


All three phases of this study further confirm through hardware 
simulation analysis that a system using a flight-qualified hand-held sextant 
and processing the observation data using statistical filter techniques have 
a basic potential of providing acceptable guidance and navigation for the 
midcourse phase. Position estimation was consistently good for all phases of 
this study, and the resulting velocity correction parameters were computed 
quite accurately from the position estimation data. The resulting capability 
of the system to guide the spacecraft to an acceptable reentry position and 
time was in all phases consistent with the requirements for a safe reentry. 

No anomalies or discontinuities were detected in the statistical filter 
processing of the hand-held sextant observation data. Operational results of 
using symbolic-alphabetic abbreviation techniques for computer control and 
display were very satisfactory. 

The use of a hand-held sextant should result in significant system 
reliability advantages. Problems due to additional error sources, such as the 
effects of sighting through a window in the space environment, variations of 
the shadows on the moon, etc., will have to be considered before actual imple- 
mentation of such a system in a spacecraft. However, the results of this 
study illustrate the basic capability of such a system to provide acceptable 
guidance and navigation information. Statistical filter theory processing of 
hand-held sextant data looks potentially good enough to warrant continued 
consideration for use in future manned spacecraft. The resulting hardware 
system configuration would offer very important advantages in simplicity and 
hence reliability. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., 94035, July 7, 1969 
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APPENDIX A 


NONLINEAR EQUATIONS OF MOTION FOR TRAJECTORY CALCULATIONS 


The equations of motion are derived under the following assumptions: 


(North) 

Z 



Sketch (a) 


1. A restricted four-body system 
is sufficiently accurate. 

2. The second harmonic term of the 
earth's gravitational field is 
sufficient . 

3. The sun and moon are spherical 
and homogeneous . 

The coordinate system is Cartesian 
and geocentric. The Z axis lies along 
the earth's polar axis, positive to the 
north. The X and Y axes lie in the 
equatorial plane with the positive X 
axis in the direction of the vernal 
equinox. The Y axis is oriented so as 
to form the right-handed orthogonal 
system shown in sketch (a) . 


The equations of motion, derived by methods given in reference 14, are: 
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where 


Re = 

/X 2 + Y 2 + Z 2 

Rm ~ 

/X m 2 + Y m 2 + Z m 2 


R s = /X s 2 + Y s 2 + Z s 2 


Am = / (X - X m ) 2 + (Y - Y m ) 2 + (Z - Z m ) 2 
As = /(X - X s ) 2 + (Y - Ys) 2 + (Z - Z s ) 2 
y = 3.986135xl0 14 m 3 /sec 2 
y m = 4.89820xl0 12 m 3 /sec 2 
y s = 1.3253xl0 20 m 3 /sec 2 

a = equatorial radius of the earth = 6.37826xl0 b m 
J = 1.6246xl0" 3 

The y values are computed by multiplying the universal gravity constant 
times the mass of the body, and J is the coefficient of the second surface 
harmonic of the earth's potential field. 

The terms involving only distances between the earth and moon, or earth 
and sun, such as y m X m /R m 3 , account for the accelerations of the coordinate 
system with respect to inertial space. 

The position of the sun and moon are obtained by interpolation of data 
from magnetic tape ephemerides. Within the sphere of influence of the moon, a 
lunar radius of 66,000 km, the origin of coordinates is translated to the 
center of the moon. Since no rotation is performed, the definitions of 
perturbations from the reference (see appendix B) remain the same. 
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APPENDIX B 


CALCULATION OF TRANSITION MATRICES 


The transition matrices used in the navigation system are obtained by 
solving linear differential equations that represent perturbations of the 
actual trajectory from the reference. Derivation of these perturbation 
differential equations is outlined below. 

The nonlinear equations of motion given in appendix A can be written in 
the form: 


X = Fi(X,Y,Z,t) ' 
Y = F 2 (X,Y,Z,t) - 
Z = F 3 (X,Y,Z,t) , 


(Bl) 


It is desired to find linear differential equations for small deviations from 
the reference. These equations may be found by expanding equations (Bl) about 
the reference trajectory in a Taylor "series and dropping all terms except the 
first order. 
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(B2) 


It is convenient to deal with systems of linear differential equations 
in multiple variables in matrix form. For this purpose, it is generally 
desirable to reduce the system to a set of first-order equations as follows. 
Define 


X 

i — 1 

n 

O) 

X 

x 4 = 6X 

x 2 = 6Y 

x 5 = 6Y 

x 3 = 6Z 

Xg = 6Z 


(B3) 
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Then the system of perturbation equations can be written in matrix form as 


dx 

dt 


F (t)x 


(B4) 


where F is a 6 x 6 matrix of coefficients and x is a 6 x 1 column vector 
of the X-j defined above. From equations (B2) and (B3) , equation (B4) can 
be written as 
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(B5) 


Consider any system of homogeneous linear first-order differential 
equations written in matrix form 


dx 

dt 


3(t)x 


(B6) 


where 3 is an n X n matrix of time variant coefficients and x is an 
n x l column vector of dependent variables. It is shown in reference 13 that 
if U is a nonsingular matrix having n columns of n linearly independent 
solutions of (B6) , then U (defined as a fundamental matrix) is a solution of 


= SCt)U (B7) 


where U(t Q ) is a constant matrix. As a special case, $ is defined as being 
the U obtained when U(t 0 ) is the unit matrix. Thus $ can be obtained 
one column at a time if equation (B6) is solved n times, each with a 
different member of x(t D ) set equal to unity and all the other members set 
equal to zero. Once $ is obtained the solution of x(t) for any given set 
of initial conditions x Q is given by 
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x(t) - $(t)x Q 


(B8) 


The matrix 0(t) represents the "transition in the states of the system of 
equations between time t 0 and t and may be written as $(t,t 0 ) to indicate 
this fact. If equation (B5) is solved in this manner, the resulting transi- 
tion matrix relates deviations from the reference trajectory at time t to 
the initial deviations at time t Q . The transition matrix $(t 2 ,ti) between 
any two times on the reference trajectory may be calculated in the same 
manner as $(t,t 0 ). 

The calculation is performed on the SDS 920 computer by solving six sets 
of perturbation equations, each with a unit initial condition on one of the 
x p, between succeeding observations. After each observation, the initial 
conditions are reset, to unity or zero, and the computation is carried out 
until the next observation. Additional details of this development may be 
obtained from reference 1. 
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APPENDIX C 


GUIDANCE, NAVIGATION, AND STATISTICAL DATA COMPUTATIONS 

LINEAR PREDICTION 


The linear prediction used is covered in detail in reference 1. It 
assumes that small deviations of the state at time t k can be obtained from 
a linear combination of the deviations at time tj c _ 1 . In particular, when 
dealing with the deviations xft]^) and x(tj c _ 1 ), the relationships may be 
written in matrix form as follows: 

x(t k ) = <Ht k ,t k _ 1 )x(t k _ 1 ) (Cl) 

x(t k ) = $(t k ,t k _ 1 )x(t k _ 1 ) (C2) 

where the transition matrix from time t k __ 1 to time t k . 

Substituting (Cl) into the definition of the covariance matrix of estimation 
error P(t k ) where 

P(t k ) = E[x(t k )x T (t k )] (C3) 

yields the following equation relating P(t k _ 1 ) to P(t k ): 

P(t k ) = Kt k ,t k _ 1 )P(t k _ 1 )$ T (t k ,t k _ 1 ) (C4) 

Similarly, substituting (C2) into the definition of the covariance matrix of 
the deviations between the actual and the reference states R(t k ) where 

R(t k ) = E[x(t k )x T (t R )] (C5) 

yields the following equation relating R(t k _ 1 ) to R(t k ) : 

R(t k ) = $(t k ,t k _ 1 )R(t k _ 1 )$ T (t k ,t k _ 1 ) (C6) 

Thus, the updating operation consists of computing equations (C4) and (C6) . 


H MATRIX 


The H matrix is computed from the reference state values in the form 


H 



8y 8y 8y\ 
3X 3Y Jz) 


(C7) 
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where the partial derivatives of the sextant angle y are found from an 
expression for y in terms of position components (X, Y, Z) and the unit 
vector to the selected star. 


KALMAN WEIGHTING MATRIX 

The Kalman weighting matrix K at time t k is computed for this 
application in the form given in reference 1, as follows: 

K(t k ) = P(t k )H T (t k )[H(t k )P(t k )H T (t k ) + Q]" 1 CCS) 

where Q is a covariance matrix of the observational errors. For this 
particular case Q is given as 


Q = [a 2 ] (C9) 

where a 2 is the variance of the error in the instrument in seconds of arc 
assuming zero mean error. 


COVARIANCE MATRIX OF THE ESTIMATION ERROR P 


The value of the covariance matrix of estimation error after a sextant 
observation P'(t k ) is computed as follows: 

P'(t k ) = P(t k ) - K(t k )H(t k )P(t k ) ' (CIO) 

This computation reflects the change in the covariance matrix of estimation 
errors P due to the observation just processed, and P'(t k ) now becomes the 
updated P(t k ) and P(t k ) is then stored. 


THE END-POINT STATE PREDICTION MATRIX A 


A prediction matrix A(t e ,t k ), which predicts the end-point state from 
the present state, is computed from 

A(t e ,t k ) = A(t e ,t k _ 1 )$(t k ,t k _ 1 )~ 1 CC11) 

An initial prediction matrix A(t e ,t Q ) is required to start the above computa- 
tation and is determined prior to the initialization of the program. To 
obtain this matrix, the system was commanded, prior to any actual data runs, 
to run from the reference trajectory injection point at t Q to the end-point 
t e without observations or velocity corrections. This resulted in a transi- 
tion matrix $(t e ,t 0 ), which was also the desired initial prediction matrix 
A(t e ,t 0 ). Updating of the prediction matrix at each observation according to 
equation (Cll) required finding $ (t k , t k _ j ) “ 1 from ®(t k ,t k _ 1 ). Because 
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*ct t j) is a symplectic matrix, it is invertable merely by rearranging 
terms and changing signs. Thus, if $(t] c ,t] c _ 1 ) is partitioned into 3x3 
submatrices 


then 








i-l 



CC12) 


(C13) 


INDICATED VELOCITY CORRECTION AV G 


With the difference between the estimated and reference states given by 
x and the updated prediction matrix from equation (Cll) the indicated veloc- 
ity correction AVq to be gained is computed from 


where G is given by 
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(C14) 


(C15) 


and Ai (t e ,t k ) and A 2 (t e ,t k ) are 3 x 3 submatrices resulting from the 
partitioning of A(t e ,t k ). 

The magnitude of AVq is 

|av g | = (av g t av g ) 1/2 


(C16) 
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COVARIANCE MATRIX AV 


The statistics of the indicated velocity correction are determined from 
a covariance matrix AV, which is computed from 

AV = E[AV g AV g T ] 

= GE [xx T ]G T (C17) 

In reference 1, it is shown that 

E[xx t ] = R - P (C18) 

so that the actual computation is 

AV = G [R - P]G t (C19) 

COVARIANCE MATRIX R 

The covariance matrix R of the errors between the actual and the 
estimated states must be corrected to reflect the effects of the velocity 
correction maneuver: The required computation, developed in detail in 

reference 1, is as follows: 

R' = (I + G) (R - P) (I + _G) T + P 
where R' is the updated R matrix and replaces R. 
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TABLE I.- INITIAL CONDITION ERRORS AT THE TRANSEARTH INJECTION POINT 

[Perilune: 70.68 hr] 


Parameter 

Initial condition errors 

X 

-0.588 km 

Y 

-.304 km 

Z 

1.254 km 

X 

.001359 km/sec 

Y 

.000949 km/sec 

Z 

.000449 km/sec 
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Mission time 
hr 


Mission function 


Observed body 


70.68 

71.0 

71.5 

72.0 

72.5 

73.0 

73.5 

74.0 

74.5 

75.0 

75.5 

76.0 

76.5 

77.0 

77.5 

78.0 

78.5 

79.0 

79.5 

80.0 

80.5 
81.0 

109.0 

110.0 
111.0 
112.0 

113.0 

114.0 

115.0 

116.0 

123.0 

124.0 

125.0 

126.0 

127.0 

128.0 

129.0 

130.0 

131.0 

132.0 

133.0 

134.0 

135.0 

144.874633 


Transearth injection time 


Observation 


Moon 


f 

Barth 


1st return velocity correction 
Observation I Moon 

Moon 
Earth 


2nd return velocity correction 


Observation 


Moon 


Earth 


( Moon 

i 

Final return velocity correction 
Time of reference perigee 
(aim point) 







TABLE III.- INITIAL CONDITION ERRORS AT 122.0 HOURS 
AS USED IN PHASE I 


Parameter 

Error 

X 

72.594 km 

Y 

1.855 km 

Z 

-19.889 km 

X 

-0.000648 km/sec 

Y 

.000328 km/sec 

Z 

.000276 km/sec 




TABLE IV.- VELOCITY CORRECTION MAGNITUDES, m/sec 



Phase I 
data 

Phase I 
Av rms (1) 

— 

Phase I 
Av rms (10) 

Phase I 
Av rms (50) 

Phase II 
data 

Phase III 
data 

Phase II 
and III 
Av rms (1) 

Phase II 
and III 
Av rms (10) 

Phase II 
and III 
Av rms( 5 °) 

First return 





2.87 

3.08 

3.35 

3.34 

3.18 

velocity correction 

— 

— 

— 

— 

±.97 

±.24 




Second return 





1.26 

1.56 

.11 

.86 

2.42 

velocity correction 

— 

— 

— 

— 

±.37 

±.76 




Third return 

0.49 

1.07 

0.81 

0.27 

.37 

.62 

.10 

.81 

2.58 

velocity correction 

±.23 




±.27 

±.30 




Total velocity 





4.19 

5.28 

3.56 

5.01 

8.18 

correction 

— 

— 


— 

±.53 



±.95 





NOTE: Plus and minus data are standard deviation data. 


TABLE V.- POSITIONAL ERROR (ACTUAL-REFERENCE) AT ACTUAL PERIGEE, km 












Maximum 

Parameter 

Phase I 
data 

^rms (1) 

*rms (10) 

^rms (50) 

Phase II 
data 

Phase III 
data 

*rms ( 1 ) 

*rms (10) 

^rms (50) 

allowable 

reentry 

corridor 

Altitude 

1.9 

±0.3 

±2.1 

±3.0 

1.9 

2.8 

±0.2 

±2.3 

±5.8 

±21.0 


±1.1 




±1.1 

±1.1 





Cross range 

.3 

±.l 

±.3 

±.4 

.2 

.4 

±0 

±.3 

±•7 

±355.0 


±.l 




±.l 

±.2 





Downrange 

-38.7 

±11.1 

±77.0 

±113 .6 

10.3 

8.9 

±8.4 

±76 .9 

±280 .8 

±806.0 


±38.8 




±38.9 

±50.9 






Single signed data are 

mean of the data 

and the plus and minus data 

are 



standard deviation data for the 

8 missions that comprised 

each data phase. 
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Figure 4.- Simulation configuration. 
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Actual velocity correction 



Figure 6.- Block diagram of digital computations. 
























Figure 7.- 


Star simulator. 


A-40576 




Figure 8.- Moon simulator 


A-40577 



A-40579 


Simulated celestial scene as viewed from the observer's station in the cab 
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Figure 10.- Lunar midcourse spacecraft cab simulator. 




Figure 11.- Computer control and display panel installed in the cab. 
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Figure 12.- Moon photograph. 
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135.0 hr 


Figure 13.- Midcourse simulation return trajectory. 




Figure 14.- Trajectory estimation errors. Phase I 
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Figure 16.- 







Magnitude of actual trajectory position errors from reference trajectory, km 


□ 261 km at 
^ perigee 



Mission time, hr 

Figure 18.- Actual trajectory position errors. Phase II. 



Magnitude of actual trajectory position errors from reference trajectory, km 
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Figure 19.- Actual trajectory position errors. Phase III. 
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